Classical simulation of entanglement swapping with bounded communication 
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Entanglement appears under two different forms in quantum theory, namely as a property of 
states of joint systems and as a property of measurement eigenstates in joint measurements. By 
combining these two aspects of entanglement, it is possible to generate nonlocality between particles 
that never interacted, using the protocol of entanglement swapping. We show that even in the more 
constraining bilocal scenario where distant sources of particles are assumed to be independent, i.e. to 
share no prior randomness, this process can be simulated classically with bounded communication, 
using only 9 bits in total. Our result thus provides an upper bound on the nonlocality of the process 
of entanglement swapping. 



By performing suitably chosen local measurements on 
an entangled quantum state, distant observers can estab- 
lish nonlocal correlations, as witnessed by the violation of 
a Bell inequality |l| . This means that quantum statistics 
cannot be simulated by classically correlated systems, un- 
less some classical communication is added to the model. 
Although experiments give strong evidence that nature 
does not use classical communication to establish correla- 
tions y , it is nevertheless interesting from a fundamental 
perspective to ask how much communication is required 
to reproduce quantum correlations. Generally referred 
to as classical simulation of entanglement, this provides 
a natural approach to the problem of quantifying quan- 
tum nonlocality. 

Nonlocality is a fundamental aspect of quantum me- 
chanics, hence quantifying it is much desirable. Besides 
being one of the most striking and counterintuitive fea- 
ture of the theory, it is also a powerful resource, allowing 
for instance for the reduction of communication complex- 
ity [3| , as well as for information processing in the 'device- 
independent' setting [^-|8[, where one wants to achieve an 
information task and prove its security without any as- 
sumption on the devices used in the protocol. 

Several works [9l4l2| underwent the task of estimating 
how much communication is needed to simulate the cor- 
relations of a maximally entangled state of two qubits 
under all possible projective measurements. This re- 
search culminated in 2003, when Toner and Bacon [13[ 
showed that one bit of comnrunication is enough. Im- 
portantly this single bit of communication is irot an av- 
erage value, but represents the exact amount that is to 
be used at each round. Thus the model is said to have 



bounded comnrunication. The communication costs of 
other states have been explored as well IJ, [la] . Notably, 
Regev and Toner [l6| have shown that the correlations 
obtainable from dichotomic measurements on any bipar- 
tite entangled state can be simulated with only two bits 
of communication, which are proven to be necessary 17 1. 
Note however that their protocol does not reproduce the 
correct marginal distributions, and that simulating more 
general measurements is much more costly in terms of 
communication [3|, [lOj. The simulation of multipartite 
entanglement also attracted some attention 18l421|. and 
two of the authors [24I recently showed that the corre- 
lations of equatorial measurements on a tripartite GHZ 
state can be simulated with 3 bits of communication, 
thus reproducing in particular the Mermin-GHZ para- 
dox, which is arguably the strongest demonstration of 
the nonlocality of this state 23|, |24 |. 

Quantum mechanics allows not only for entangled 
states of distant systems, but also for entangled mea- 
surements. In such a measurement the initial state is 
arbitrary — it could be entangled or not — but the final 
state is entangled, that is the eigenstates of the operator 
that represents such a measurement are entangled. This 
second aspect of entanglement is in itself independent of 
nonlocality — although it leads to nonlocality when com- 
bined with entangled states [25| . It demonstrates another 
nonclassical feature of entanglement, which is, loosely 
speaking, the possibility to ask two (or more) quantum 
systems questions about their relations without gaining 
any information about the individual properties of each 
subsystem [26j . 

It is therefore a natural question to ask whether or 



not the classical simulation of protocols involving both 
entangled states and entangled measurements is possible 
with bounded communication, and how much communi- 
cation is required. Here, we investigate this question for 
the scenario of entanglement swapping [27| , where quan- 
tum particles that never interacted become nonlocally 
correlated after their twins underwent a joint measure- 
ment. We work in the scenario of bilocality (28|, |29|, in 
which the shared randomness we will use in the simula- 
tion protocol is assumed to originate from two indepen- 
dent sources (see figure [T]). Since entanglement swapping 
can be achieved with fully uncorrelated quantum sources, 
even experimentally [30, |3l| , it is indeed natural to im- 
pose an equivalent constraint on the simulation model. 

Importantly the constraint of bilocality makes the sim- 
ulation of entanglement swapping a challenging problem. 
The goal is basically to generate singlet nonlocal corre- 
lations between two parties which are initially fully un- 
correlated. At first sight this may even appear to be 
impossible with finite communication, in the light of a re- 
sult by Massar et al. [32| , showing that simulating singlet 
correlations requires either infinite shared randomness or 
infinite communication. Here however, we will show that 
the process of entanglement swapping can be simulated 
with 9 bits, by presenting an explicit protocol. Since 
the communication cost of the cheapest protocol can be 
considered as a measure of its nonlocality, our result pro- 
vides an upper bound on the nonlocality of entanglement 
swapping. 

The entanglement swapping process. We consider three 
distant parties, Alice (A), Bob (B) and a Referee (i?). R 
shares two maximally entangled qubit pairs in the state 
|*_) = -75(101) - |10)) with A and B, respectively. The 
first pair is produced by a source located between A and 
R, the second one is produced by an independent source 
located between B and R, see Fig. [1] (i). Accordingly, 
Alice and Bob are initially uncorrelated. By performing 
a Bell state measurement, i.e. a two-qubit joint measure- 
ment which features four maximally entangled Bell states 
as cigenstates, R projects A and B's particles onto one of 
the Bell states. R can thus 'swap' entanglement to A and 
B. This protocol is essentially identical to the celebrated 
quantum tclcportation protocol [33| ; entanglement swap- 
ping is basically a teleportation of entanglement. 

To complete the protocol, the Referee needs to com- 
municate the result of his measurement r — by sending 
two bits of classical communication — to (say) Bob, who 
can then apply a suitable unitary local transformation to 
his qubit to finally share a definite Bell state with Alice, 
say 1^1'-). Upon receiving measurement settings x and 
y (represented by unit vectors on the Bloch sphere) and 
performing the corresponding projective measurements 
on their respective qubit, Alice and Bob obtain binary 
measurement outcomes a = ±1 and 6 = ±1 respectively. 



Alice 



Referee 



Bob 



(i) 



X 

\ 

n 

\ 

a 



(ii) X 



\'^~)ar 
— * 



|*-)i?B 
* 



i 



— * — 



■ * — 



y 

I 

D 

I 

b 

y 

i 

n 
b 



FIG. 1. (i) The scenario of entanglement swapping with two 
fully independent sources of |^-) states. The Referee sends 
to Bob the result r of his Bell state measurement. Upon 
receiving these two bits of communication, Bob can apply the 
adequate unitary operation to his qubit such that Alice and 
Bob finally share a singlet state, (ii) The classical simulation 
of entanglement swapping in the bilocal scenario, i.e. with 
two uncorrelated sources of shared random variables Xar and 
Xrb ■ The three parties exchange messages symbolized by the 
thick dashed arrows. Here we present a simulation protocol 
using 9 bits of communication in total. 



These outcomes exhibit nonlocal correlations of the form 

E{x^ y) = P{a = b\x, y) — P{a ^ b\x, y) = —x ■ y, (1) 

with random marginals. 

A classical simulation of this protocol in the bilocal sce- 
nario would then amount to the following (see Fig.[l](ii)): 
Alice and Bob receive as inputs the measurement direc- 
tions X and y, independently of the random variables 
Xar and Xub they may each share (independently) with 
R. After a finite amount of information exchange be- 
tween the three parties, Alice and Bob produce their 
outputs which, after averaging over the shared random 
variables, must be correlated as in ([T]) j34| . 

At this point, it is instructive to recall the result of 
Massar et al. 34] , which implies that correlations of the 
form ([T]) are impossible to simulate with finite commu- 
nication when A and B only have finite shared random- 
ness (in the sense that their shared randomness could be 
established with bounded communication) — even if they 
may each have an additional source of infinite random- 
ness, independent for A and B. The present scenario 
may at first glance appear similar: initially, A and B 
each only have access to independent sources of infinite 
randomness, Xar^ and Xrbi respectively. After any fi- 
nite exchange of communication between the parties (in- 
cluding R), they could only share some finite amount of 
randomness. Intuition suggests that the result of [32[ 
implies that, from such randomness, correlations of the 
form II]) cannot be simulated with finite communication, 
and entanglement swapping is impossible to simulate in 
a bilocal manner 13 5| . However, the present scenario dif- 
fers from that of [32| in that the third party, R, can 



share an infinite amount of randomness with A and B 
independently — by having aceess to both Xar and Xhb — 
and can exchange some (finite number of) bits with A and 
B. This seemingly subtle difference dramatically changes 
the situation: despite the above heuristic reasoning, en- 
tanglement swapping can be simulated with finite com- 
munication in a bifocal manner, as we now show. 

Simulation protocol. We start by deriving a proto- 
col for the simulation of equatorial measurements x = 
(cos(/)A,sin(/)A,0) and y = (cos 0b, sin (/j^, 0). In this 
case, the correlation ([T|) is equivalent to 



P(a= 



1— cosf 



»Bj = 



= sm 



9s \ 2 



(2) 



Note in particular that if (f>A ~ (I^b, then P{a ~ b\(j)A = 
(J)b) = 0, as expected from the anticorrclations of the 
singlet state. 

Let us divide the equator of the Bloch sphere into 2m 
sectors of the same size, and assume that Alice and the 
Referee (resp. the Referee and Bob) share a random vari- 
able Xar. (resp. Xrb)i with Xar. and X^b uniformly dis- 
tributed in [0, —1. Let Alice tell Bob in which sector her 
setting (modulo tt) lies, and send him the additional bit 
ca = [{4>A rnod — ) < Xar\ (we use the Iverson bracket, 
with [E] = 1 if the expression £ is true, [£] = otherwise); 
let the Referee send to Bob the bit cr = [Xar. < Xrb]; 
and let Bob calculate cb = [Xrb < {4'b mod ^)]. 

To warm up and gain some intuition, assume that Alice 
and Bob's settings both lie in [0, — ], and let Alice simply 
output a = +1. li Ca = Cr = Cb, then Bob knows 
that (with probability 1) (f>A ^/^ 4>b', let him then also 
(for now) output b — a — +1. On the other hand, if 
Ca 7^ Cr or Cr ^ Cb, it could be the case that (f>A ~ 4'b', 
to ensure that P{a = b[(t>A = 4'b) = 0, Bob must output 
6 — —a = — 1. The probability that Alice and Bob give 
the same output is then (assuming 4a < 4b', the case 
4'A > 4'B being similar) 
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If "1 > ^! this probability is too large (for any 
4'A, 4'B & [0' ml)' compared to the desired value of 
^"^( 2 ) ®' -Now, one can actually make use of 
some freedom in the case ca = cr = cb, and ask Bob to 
output b = a = +1 only with some probability p (and to 
output 6 = — a = — 1 with probability 1 — p). By letting 
p depend on 4>b and Xrb, and take the particular form 
p{4>b-Xrb) = 2^ cos{4)b->^rb), we obtain, as desired. 
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In order now to extend this protocol to all possible 
equatorial measurements 4'A, 4'b S [0,27r], one needs to 
find adequate strategies for Bob (i.e. adequate functions 
p{4'b—Xrb) G [0, 1]) for all possible sectors where Alice 
and Bob's settings may lie, and for all possible values of 
ca , Cr. and cb ■ We were able to find such strategies for 
the case 777 = 4, see Table ID This leads us to define the 
following protocol: 

Protocol 1: Let Alice and the Referee (resp. the Ref- 
eree and Bob) share a random variable Xar (resp. Xrb), 
with Xar and Xrb uniformly distributed in [0,j]. Af- 
ter reception of Alice and Bob's measurement settings 
4'A, 4'B, the three parties proceed as follows: 

• Alice calculates a = sign(sin0^) = ±1, ja ~ [-{4'A 
mod7r)J e {0,1,2,3}, 0^ = {4)a-Jaj mod tt) e 
[0,j[ and CA ~ [4>'a < ^ar] S {0,1}. She sends ja 
(2 bits) and ca (1 bit) to Bob, and outputs a. 

• The Referee calculates cr = [Xar. < Xrb] G {0, 1} and 
sends it to Bob (1 bit). 

• After reception of ja, Bob calculates (3 = sign[sm.{4>B ~ 



JA^)] 



b'B 



f>B - JAj modTr) e [0,7r[. He 



determines the index jb = 1^4''b\ ^ {0,1,2,3} 
of the sector where his angle 4''b ^'^^^' '^'"^^ the bit 
Cb ~ [^RB < 4^'b^Jbj] G {0,1}. Depending on 
JB, CA, Cr, cb,4''b '^'^'^ ^RB , he outputs h = j3 with prob- 
ability pi'^CRCB('^'B~^RB) , and b = —/3 with probability 
1 - PilcnCB^^B-^^RB), for the functions pi%j^^^ de- 
fined in Table\^ 



As explicitly shown in the Appendix below, the above 
protocol gives the desired probability P(a = &|0a, 4'b) = 
sin(-^ii^^)2 for all possible equatorial measurements, us- 
ing 4 bits of communication. It can then be extended in 
the following way to all measurement directions on the 
Bloch sphere, using a similar technique as in [lO| : 

Protocol 2: For measurement directions 

X ~ (sin6'ACOS0yi, sin^A sin0A,cos0A) and y = 
(sin 9b COS 4'b , sin 9b sin 4>b , cos 9b ) , the three parties 

• run Protocol 1 with input angles 4>A and 4>b', Alice and 
Bob obtain intermediate outputs qq and bg. 

• run Protocol 1 a second time (using a new set of in- 
dependent variables Xar, Xrb), now with input angles 
aQ9A and —b^OB; Alice and Bob output the outcomes a 
and b of this second run of Protocol 1. 

This second protocol now simulates the desired correla- 
tion E{x, y) = X ■ y for all possible projective measure- 
ments by Alice and Bob, with 8 bits of communication; 
for more details on the calculations, see the Appendix. 
Note that Protocols 1 and 2 do not simulate the correct 
marginals. In order to randomize the marginals, Alice 
can — at the very end of the protocol — generate a ran- 
dom bit and send it to Bob; depending on the value of 
this bit, they will both flip their outcomes or not. All in 
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^(7), for all values of Jb,ca,ch and cb, defining the probability p\ 
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outputs 6 = /3 in Protocol 1. Note that in all cases, 
interval. 



,(7) G [0, 1] for all possible values of 7 = (J)'b — ^rb in the given 



all, the entanglement swapping correlations can thus be 
simulated with 9 bits of communication. 

Discussion. We thus have proved that remarkably, 
the entanglement swapping process can be simulated 
with bounded communication, even in a bilocal scenario 
where Alice and Bob are (as in the quantum case) com- 
pletely uncorrelated before the protocol is run, and there- 
fore do not have any prior shared randomness. Our pro- 
tocol provides an upper bound on the nonlocality of en- 
tanglement swapping in terms of its communication cost. 
It is an open question whether fewer bits of communi- 
cation are actually sufficient: it might indeed be possi- 
ble to simulate equatorial measurements more efficiently 
than with Protocol I, or to find a more direct protocol 
that docs not treat separately the azimuth and zenith 
angles of the measurement settings, more in the spirit of 
the Toner-Bacon simulation protocol for singlet correla- 
tions [13|. 

Next, it is natural to consider the simulation of multi- 
stage entanglement swapping, which is essential for long 
distance quantum communication. Now, N referees (i?i, 
i?2, ..., Rn) are placed on a line between Alice and 
Bob. Two neighboring referees share a singlet state, 
while i?i and Rn share singlet states with A and B, 
resp.; each referee performs a joint measurement, leav- 
ing at the end the particles of Alice and Bob entan- 
gled. Whereas the quantum protocol has a straightfor- 
ward and nice iterative character, we were not able to find 
a simulation protocol with a finite amount of communi- 
cation in a (iV-l-l)-locality scenario 129|. Consider for 
instance the case with one additional referee i?2- Anal- 
ogously to our Protocol 1, assume that Ahce and i?i 
share the random variable Xar^ Ri and R2 share XFI1R21 
and i?2 and Bob share Xr^b, all uniformly and indepen- 
dently distributed on some interval [0,—]. After some 
finite communication. Bob could for instance (as in our 
first attempt, before Eq. ([3])) output fe = a = I if and 
only if [(j)A < \ARt] = [AaAi < ^RirA = [AfliJ?2 < 
Xr2b\ — [XR2B < 'Pb]- This would result in the proba- 



bihty P{a = h\4>A,^B e [0, ^]) = p\^A - M^ which 
scales cubically with (pA — (f'B , and is therefore too small 
when (j)A is close to (f>B- It is unclear how to change the 
cubic scaling with finite communication . The follow- 
ing questions remain open; can multistage entanglement 
swapping be simulated with finite communication? Or 
can one prove, that above a certain value of A^, an infi- 
nite amount of communication is necessary? 
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APPENDIX 

In this appendix we prove that Protocols 1 and 2 
generate the desired correlations ([T|), or in other terms, 

P[a = b\x,y) = 



Bob's side, ensure that the following relations hold, as 
required: 



_ l-xy 

2 



Correlation from Protocol 1 



We first note that the definitions of a = sign(sin(/)^) 
and of (/)^ = {<pA—JA'^ niod tt) G [0, f [ on Ahce's side, 
and the definitions of /? = sign[sin(0B — jyii-|)] (and then 
of 6 = ±/3) and of (jy'^ ~ (i/ib — ja| mod tt) e [0, 7r[ on 



P{a = 6I0A, 4>b) = P{a ^ h\4>A, 4>b + tt) 

= P{a ^ b\(j)A + TT, 4>b) = P{a = h\(j)A + tt, 0b + tt) 

TT TT 

= P{a = b\(l)A+i-,4>B+i^) for any j e Z. 



It is therefore sufficient to check that the correct correla- 
tion is obtained for 4>a £ [0, f [ and 4>b G [0,7r[. 

For such values of 0^, 0_b (for which (j)'^ = 4>a and 
4''b ^ 4'b), the probability P{a = h\(l)A,4>B) obtained 
from Protocol 1 can be calculated as follows: 
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= — I / dXAR / dXRB Pooi{<Pb->^rb) 



jBii rrpB-jB-i 

dXAR / dXRB Pon(0B-AK_B) 

J\ar 

<t>B-JBf nf 

dXAR / dXRB Poio(<^s-A_rb) 

JlpB-JB^ 
<^A /••I^B-JB^ 

dXAR / rfA_RB Pooi(<?^b-A_rb) 

dXAR / dAfls pioai^B-XRB) 

<B-JB^ JtpB-JB^ 

<t>A l-f 

dXAR / dXRB Poio(^b-Ahb) 

4>B-JBJ J >^AR 

f<pB-JE ' 

dXAR ^\^^,^^ 

Jo 

_; \ , 7 R / i 

0A J<t>B~JB^ 



dXRB p{^i{(t>B-XRB) 



dXAR I dXRB p{oo{4>b-Xrb) 



dXAR / dXRB Pnoi(f>B-XRB) 

4>A J >^AR J 



One can then check that with the choice of functions 
^ACflCs G [0,1] indicated in Table H this leads (for all 



as desired. 



Correlation from Protocol 2 



After running Protocol 2 for inputs x 
(sin6'yi cos0A,sin6'A sin(/)yi, cos6'a) and y 
{sin 9 B cos (j) B, sva 9 B sin 4>B, cos 6 b), the probal 
that Alice and Bob's outputs are the same is 



,M 



P{a = h\x,y) 

= ^ P[ao,bii\(j)A,<j)B)P[a^h\6A,9B,ao 

ao,fco = ±l 

E„, , I , , ■.l~cos{aoBA+bt)9B) 
P[ao, bo\(j)A,(PB) ^ 

ao,fco = ±l 

= P{ao = bo\cpA,(PB) ^ 

^ p{ -^h\A ^ ^ t - cos{9a-6b) 
+ P(ao / bolcpA, 9b) 2 

_ 1 - cos{4>a-4>b) 1 - cos(6'a+6'b) 
~ 2 2 

1 + cos((/)a-(/'b) 1 - cos(6>a-6>b) 
2 2 

1 — cos 6*^ cos 6*5 — sin0Asin0BCOs(0A— 0s) 
2 
_ 1 -a; ■ j/ 
2 

Protocol 2 thus reproduces the desired entanglement 
swapping correlation ([T]). 



